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ABSTRACT 

We have investigated the problem of sunspot magnetohydrostatic equihbrium with comprehensive 
IR sunspot magnetic field survey observations of the highly sensitive Fe I lines at 15650 A and nearby 
OH lines. We have found that some sunspots show isothermal increases in umbral magnetic field 
strength which cannot be explained by the simplified sunspot model with a single-component ideal gas 
atmosphere assumed in previous investigations. Large sunspots universally display non-linear increases 
in magnetic pressure over temperature, while small sunspots and pores display linear behavior. The 
formation of molecules provides a mechanism for isothermal concentration of the umbral magnetic 
field, and we propose that this may explain the observed rapid increase in umbral magnetic field 
strength relative to temperature. Existing multi-component sunspot atmospheric models predict that 
a significant amount of molecular hydrogen (H2) exists in the sunspot umbra. The formation of H2 can 
significantly alter the thermodynamic properties of the sunspot atmosphere and may play a significant 
role in sunspot evolution. In addition to the survey observations, we have performed detailed chemical 
equilibrium calculations with full consideration of radiative transfer effects to establish OH as a proxy 
for H2, and demonstrate that a significant population of H2 exists in the coolest regions of large 
sunspots. 

Subject headings: magnetic fields — molecular processes — Sun: infrared — Sun: sunspots 
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1. INTRODUCTION 

1.1. Background 

Sunspot umbrae are quasi-stable structures that can 
be considered to be in nearly magnetohydrostatic (MHS) 
equilibrium due to their long lifetimes relative to the dy - 
namical timescale of the quiet sun (Meyer et al. 1977). 
Hale's historic observation of the Zeeman effect in 
sunspots revealed them as concentrations of strong mag- 
netic fields (jHald [ToOSh . and Biermann and Alfven pro- 
vided us with the first phy sically plausible explanation o f 
the sunspot phenomenon (Biermannl l 19411 : lAlfven|[l943[ ). 
Although the sunspot atmosphere is cool and has a 
low degree of ionization, its conductivity is high, and 
it can be assumed that the magnetic field is "frozen 
in" to the gas. In Alfven's theory the strong magnetic 
field suppresses the convective heating of the interior of 
the sunspot. It also supports the cool sunspot interior 
against the higher pressure of the hotter surrounding 
quiet sun (Deinzer 1965). For a sunspot with circular 
symmetry the MHS equilibrium state at any radius from 
the center of the sunspot (r) and height (z) in the sunspot 
atmosphere can be written as: 



%{z)=n,{r,z)m{r,z) + ^{BUr,z)+F,ir,z)), (1) 



^ Visiting Student, National Solar Observatory, operated by 
the Association of Universities for Research in Astronomy, Inc. 
(AURA), under cooperative agreement with the National Science 
Foundation. 



for c.g.s. units, where Pqs is the gas pressure in the quiet 
sun, and the pressure in the sunspot is given by the prod- 
uct of the number density of the gas (n^), Boltzmann's 
constant (A:), and temperature (Tg). The magnetic field 
contributes a pressure term from the vertical component 
of the field {Bz) and a tension or curvature term (Fc) 
which is given by: 



Fc{r,z) 



'^.(r^.)^^^^r^ (2) 



where the integral runs from the radius in question 
to the maximum radius of the sunspot (a), and re- 
quires knowledge of the vertical gradien t of the radial 
compo nent of the magnet ic field (Br) (iCowlingI 119761 : 
Marti nez Fillet Vazquezlll993f ). 

Although this theory provides us with a basic descrip- 
tion of the MHS equilibrium condition for sunspots, ob- 
servational tests are difficult to carry out, even for the 
form of circular symmetry assumed in Equation [TJ Be- 
cause of radiative transfer (RT) effects, the magnetic and 
thermodynamic properties of a sunspot observed in re- 
gions of different temperature do not originate at the 
same geometrical depth in the sunspot atmosphere, and 
information about the vertical structure cannot be de- 
rived easily. Therefore the contribution of the curvature 
force term, which requires knowledge of the vertical gra- 
dient of the magnetic field, cannot be determined directly 
from observations. Nevertheless, under the restrictive 
conditions present in the sunspot umbra it is possible to 
make some simplifying assumptions. Here the magnetic 
fields are mostly vertical and the contribution from the 
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Fig. 1. The B'^ vs T relation from NOAA 11130 observed with 
FIRS on December 2, 2010 displaying an isothermal increase in 
the umbral magnetic field. Points selected from the penumbra are 
shown in gold (+) and points from the umbra are black (x). 



curvature force term is expected to be suiall. In addi- 
tion, the radial gradient of temperature is small, and the 
magnetic field and temperature observations can be con- 
sidered to originate at a single height. Therefore, if we 
make the assumptions that sunspots are vertical mag- 
netic fiux tubes with constant magnetic field strength, 
and have an ideal gas atmosphere with constant den- 
sity and height, then Equation [T] can be simplified to the 
thermal-magnetic relation: 

B\r)^T,,-T,{r). (3) 

Although Equation [3] is elegantly simple it has not 
been verified observationally. While some examples of 
the empirical relation between and T determined 
from sunspots are best described by a linear relationship 
(iGurman fc House! 119811 : iLivingstonI I2QQ2I : iPenn et alJ 
I2QQ2L l2QQ3af ). the majority of modern observations 
show that can be a hig h ly no n- hnear function 
of T (Balthasar fc Schmidt! Il993l; ! Kodd fc Rabin 
19931 : iLites et alJ 11993!: !Martmez Pill et fc Vazquez 
1993'; 'Solanki et alJ !l993!: iStanchfield^et alJ !l997r 
Westendorp Plaza et aD !2QQl! : !Penn et alJ !2QQ3b! : 
Mathew et al. 2004). Figure [T] gives an example of 
the thermal-magnetic relation for the main sunspot 
in NOAA 11130 that we have recently derived from 
spectropolarimetric observations at 15650 A, where 
temperature has been determined from the contin- 
uum contrast, and magnetic field was obtained by 
Milne-Eddington inversion of the Fe I line pair. The 
non-linearities in the B^ vs. T relation are evident: the 
sunspot shows an unexpected sharp increase in magnetic 
field strength at relatively constant temperature in the 
darkest part of the umbra. This result is puzzling. 
With the simplifying assumptions leading to Equation 
m Figure [T] implies the existence of a heretofore uniden- 
tified mechanism at work in the coolest regions that 
reduces the thermal pressure of the sunspot atmosphere 
to accommodate the higher magnetic pressure in these 
regions. 



1.2. The Magnetohydro static Equilibrium of Sunspots 
and the Role of Molecular Hydrogen 

In deriving Equation l3] we have made the assumptions 
that: 1) sunspot magnetic fields are near vertical with 
respect to the local solar surface and that the curva- 
ture force has a negligible effect on the force balance of 
a sunspot, 2) the geometrical depth of the r = 1 level 
depends roughly on the temperature T, and that the 
temperature and magnetic field observed in a sunspot 
originates from roughly the same height in the sunspot 
atmosphere, and 3) the density of the solar atmosphere 
is a constant that does not change with temperature. In 
assuming a fiux tube model, we also neglect fine struc- 
ture, flows, growth, decay, and other evolutionary effects. 
While the first two assumptions are diflicult to verify due 
to radiative transfer effects, they may be justified under 
restrictive conditions such as those encountered in the 
sunspot umbra. On the other hand, the third assump- 
tion of constant density is obviously not warranted, espe- 
cially under the temperature regime encountered in the 
umbral photosphere. 

At an effective temperature of 5,750 K the quiet sun 
photosphere is largely neutral and already harbors many 
heavy-elem ent molecular species w ith high dissociation 
energies (Grev esse fc Sauva ililll. At even lower tem- 
peratures in the sunspot umbra many more molecular 
species, in particular molecular hydrogen (Il2), are able 
to form. Because the abundances of heavy elements in 
the solar photosphere are very low compared to hydro- 
gen, the number of molecules formed containing heavy 
elements is correspondingly low with respect to the to- 
tal particle number density, therefore the formation of 
heavy element molecules cannot have a significant ef- 
fect on the thermodynamic properties of the solar atmo- 
sphere. However, in the umbra the atmosphere is very 
cool, making it possible to form a substantial fraction of 
H2. Observations of the fiuorescent H2 lines in the ultra- 
violet have confirmed its presence in the chromosphere 
above sunspots (Jord an et all 11978! : !Bartoe et al.! !1979! : 
!lnnes! |2008>) . while atmospheric models of the sunspot 
umbra predict a molecular hydrogen population of up 
to 10%, peakin g near the height of continuum formation 
(jMaltbv et al.! [l986). 

The formation of a large fraction of molecules may have 
important effects on the thermodynamic properties of 
the solar atmosphere and the physics of sunspots. For 
example, as free atoms combine into molecules the dis- 
sociation energy is released. This may increase the local 
thermal energy content of the atmosphere if it cannot be 
dissipated rapidly enough. On the other hand, molecules 
have the ability to store energy in rotational and vibra- 
tional degrees of freedom which do not contribute to the 
thermal signature of the gas, resulting in an increased 
heat capacity. Therefore if a significant fraction of the 
hydrogen atoms in the photosphere of a sunspot exist 
in molecular form, then the thermodynamic properties 
of the sunspot will be significantly different from that 
of the quiet sun photosphere. Finally, the combination 
of free atoms into molecules decreases the total particle 
number density and pressure of the gas. This effect is 
of particular importance to the problem of MHS equilib- 
rium in sunspots since it provides a mechanism for chang- 
ing the gas pressure of the sunspot atmosphere without a 
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TABLE 1 
Summary of Selected Observations 



Target Name 


Date 


UT Time^ 


Lat [°] 


Lon [°] 




Classification 


Phase 


NOAA 9429^ 


2001-04-18 


14:44:17 


8.39 


13.42 


0.947 


OL 


decay 


NOAA 11035 


2009-12-17 


15:58:13 


27.73 


28.83 


0.754 




growth 


NOAA 11046 


2010-02-13 


17:30:53 


22.88 


3.41 


0.859 




decay 


NOAA 11049 


2010-02-19 


15:57:02 


-20.61 


17.72 


0.933 




quiescent 


NOAA 11101 


2010-09-02 


14:27:11 


10.97 


38.70 


0.782 


a 


quiescent 


NOAA 11130 


2010-12-02 


15:42:30 


11.03 


49.91 


0.624 




decay 


NOAA 11131 


2010-12-06 


17:00:35 


29.78 


-18.45 


0.815 


OL 


quiescent 



^ Start time of the observation. 
^ Cosine of the hehocentric angle. 

^ Preliminary observation taken with the Horizontal Spectrograph. 



corresponding change in the temperature, which may ex- 
plain the isothermal intensification of the magnetic field 
strength in the darkest regions of the sunspot seen in 
Figure [H Even a small alteration in the balance of gas 
pressures in a sunspot may translate into large changes 
in the sunspot magnetic field. 

1.3. Objectives of This Research 

Given the critical role that H2 may potentially play in 
the physics of sunspots, it is the objective of this study 
to determine if a significant fraction of H2 exists in the 
photosphere of sunspots. We will also demonstrate that 
H2 formation occurs at a temperature coinciding with the 
sharp rise of B'^ seen in Figure [T] to support the argument 
that this feature is due to the formation of H2 . 

1.4. Methodology 

To establish the existence of H2 in the sunspot photo- 
sphere and its link with the isothermal intensification of 
the umbral magnetic fields, it is necessary to characterize 
the 2D magnetic field, temperature, and H2 fraction of 
many sunspots at different stages of evolution. In pre- 
vious investigations, the authors presented the thermal- 
magnetic relation from a single sunspot observation, or 
from observations of many sunspots with very limited 
spatial resolution. To fully explore the relationship be- 
tween temperature and magnetic field strength, we have 
conducted a comprehensive survey of sunspots, obtaining 
high spectral and spatial resolution spectropolarimetric 
observations of the infrared Fe I line pair at 15650 A. The 
instrumentation and observational setup, as well as the 
dataset, are described in 321 The reduction of the data, 
including the treatment of scattered light and the mea- 
surement and calibration of the telescope and instrumen- 
tal polarization crosstalk, will be described in ^3.11 and 
13. 2[ We have derived the vector magnetic field configu- 
ration of the sunspots using a Milne-Eddington inversion 
code (3331). 

Due to the highly forbidden nature of the infrared ro- 
vibrational transitions of H2 in the solar atmosphere, 
and the current lack of space-based instruments to pro- 
vide spectral coverage of the fiuorescent emission lines 
in the UV, direct measurement of the abundance of 
H2 abundance in the photosphere lies beyond our grasp. 
Nevertheless, since the formation of the molecules is di- 
rectly controlled by the temperature of the gas, it is 
possible to infer the photospheric H2 abundance from 
proxy measurements of other similar molecules. The dis- 
sociation energy of the H2 molecule (4.48 eV) is very 



similar to the dissociation energy of the hydroxide (OH) 
molecule (4.39 eV) which exhibits strong absorption fea- 
tures throughout the visible and IR sunspot spectrum. 
Serendipitously, it has strong absorption lines inside the 
spectral window of the Fe I line pair at 15650 A used for 
the measurement of sunspot magnetic fields. We have 
established the equivalent width of the OH 15652 A line 
as a proxy for the H2 abundance through detailed calcu- 
lation of the chemical equilibrium, and OH spectral line 
synthesis for model solar atmospheres with different ef- 
fective temperatures. The details of this calculation are 
described in ^ 

We will present the results of this study, including 
the measurement of the magnetic field pressure 5^, the 
equivalent width of the OH 15652 A line, and the in- 
ferred H2 abundance as a function of the temperature of 
the sunspot atmosphere in ^Sl Finally, we will summa- 
rize this study and discuss the influence of H2 formation 
on the emergence and evolution of sunspots in ^ 

2. OBSERVATIONS AND INSTRUMENTATION 

The data present in this paper were obtained with 
the Dunn Solar Telescope (DST) at the National So- 
lar Observatory (NSO) located in Sunspot, New Mexico. 
Preliminary observations were taken with the Horizon- 
tal Spectrograph (HSG) in 2001 and 2005, and the main 
survey observations were taken over a period from July 
2009 to December 2010 using the newly constructed Fa- 
cility IR Spectropolarimeter (FIRS). For every sunspot 
in the survey we have made spectropolarimetric observa- 
tions of a 10 A bandpass at 15650 A covering the Fe I 
15648.5 A g=3 and Fe I 15652.9 A g=1.7 lines, a blend 
of two OH lines at 15650.8 A, and two OH lines on ei- 
ther side of the Fe I g=1.67 line at 15651.9 and 15653.7 
A. The complete survey contains 66 observations of 23 
different active regions and pores. From this sample we 
present seven illustrative cases. Table [1] summarizes the 
observed active regions, the date and time of observa- 
tion, and their heliocentric position. Th e full su r vey w ill 
be presented and discussed at length in iJaeg gl? (20111). 

The Horizontal Spectrograph (HSG) is a versatile and 
reconfigurable 3 m spectrograph arm with multiple avail- 
able gratings. It consists of an adjustable rectangular 
slit, collimator lens, plane grating, camera lens, and de- 
tector. The experimental HSG setup for thes e sunspot 
observ ation s wa s similar to those described in iLin et all 
(|1998[ ) and [03 (|1995f ). During the 2001 observations 
a 256 X 256 HgCdTe NICMOS 3 detector was used to 
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Fig. 2. — Comparison of data from Hinode SOT/SP and FIRS, using coincident sunspot observations taken on 2009-07-07. Shown are 
the intensity maps, full Stokes spectra, and line spectra from FIRS (red) and Hinode (black) at the positions indicated by the dotted lines. 



record Stokes spectra. 

FIRS is a new facility instrument for the DST recently 
completed by the Institute for Astronomy at the Univer- 
sity of Hawai'i (P.I. Haosheng Lin) in collaboration with 
the NSO. Several unique design elements allow FIRS 
to perform full- Stokes vector spectropolarimetry at two 
wavelengths simultaneously in the visible and infrared, at 
6302 and 10830 A or 6302 and 15650 A. An off-axis near- 
Littrow configuration and steeply blazed grating provide 
high spectral resolution, and spectra are imaged using a 
Kodak 2048 x 2048 CCD for the visible and a Raytheon 
Virgo 1024 x 1024 HgCdTe array for the infrared. FIRS 
maps four slit positions simultaneously to increase ob- 
serving cadence and make efficient use of the large format 
of these detectors. A scan area of 150 x 75 arcsec^ can 
be covered in 20 minutes with 0.29 arcsec^ spatial pixels, 
while achieving a signal-to-noise of 1000. Specialized nar- 
row band filters for each wavelength prevent the spectra 
from adjacent slits from overlapping. Diffraction-limited 
performance at 6302 A can be achieved with FIRS using 
the DST^s Hi^h Orde r Adaptive Optics (HOAO) system 
(|Rimmele et al.|[2004f ). which also provides excellent im- 
age correction at infrared wavelengths even in mediocre 
seeing conditions. Separate liquid crystal variable re- 
tarders for each wavelength modulate the polarization 
in an efficiency-balanced scheme which is analyzed and 
split into two orthogonally polarized beams by a Wol- 
laston prism which allows the spurious seeing effects to 
be eliminated in post-processing. A more detailed de- 
scription of the instrument can be found in i Jaeggli et al] 
(j2010.) and soon in Lin et al. (in prep.). 

3. DATA REDUCTION AND ANALYSIS 

3.1. Data Reduction 

The observations from FIRS were processed using the 
IDL reduction tools (http://kopiko.ifa.hawaii.edu/firs/) 
which corrected for the infrared detector non-linearity. 



subtracted darks, and demodulated the raw polarization 
states to retrieve the Stokes vector. A fiat field was pro- 
duced by fitting and removing absorption lines from spa- 
tially averaged quiet sun observations and was applied 
to Stokes I to remove intensity variations due to the de- 
tector, slit, narrow-band filter, and other elements in the 
optical system. The demodulation scheme for the Stokes 
vector effectively removes response variation under the 
assumption that the filter profile is fiat over the spectral 
window. The spectra from each slit were corrected for 
spectral curvature and drift, the orthogonally polarized 
beams were combined, and the spectra were stacked into 
a data cube. Following the reduction the data cubes were 
cropped spatially to a relevant area around each sunspot. 

Observations from the HSG were reduced in a simi- 
lar fashion to th e FIRS observations (details in lLin et al.l 
(1998) and iLinI ()1995f )). however no correction was ap- 
plied to compensate for the detector linearity. 

3.2. Correction of Instrumental Polarization 

The Stokes Q, U, and V spectra were corrected for 
in strumental polariz ation cross-talk using the method 
of iKuhn et al.l ()1994f ). This approach rather simply as- 
sumes that the Stokes Q and U profiles should be sym- 
metric about the line center and Stokes V should be anti- 
symmetric. From this assumption the cross-talk coeffi- 
cients are determined from a small portion of the spec- 
trum and applied to the full data cube. 

Based on these assumptions, one would think that the 
Kuhn method might not work well for physical situa- 
tions which produce non- symmetric and anti-symmetric 
Stokes profiles. Therefore, we have verified the appro- 
priateness of the Kuhn method, and the data quality in 
general, by comparing corrected FIRS observations to re- 
duced spectra from the Hinode Solar Optical Telescope 
Spectropolarimeter (SOT/SP) which were closely coin- 
cident in time (about 4 minutes apart). The polariza- 
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Fig. 3. — Point -by-point scatter plot comparison of the full Stokes spectra from Hinode SOT/SP and FIRS data cubes, using coincident 
sunspot observations taken on 2009-07-07. The color scale indicates the density of points, 50, 75, 95% contours are displayed for I and 90, 
95, 99% contours for Q, U, and V. The solid line indicates the 1:1 correspondence. 



tion cross-talk in the SOT/SP has bee n very carefully 
measured and corrected as described by llchim oto et aL 
(^2008). The maps and example spectra from each obser- 
vation are shown in Figure [21 and in Figure [3] we show a 
point-by-point comparison of all quiet sun and sunspot 
spatial and spectral positions from the aligned datacubes 
of Stokes spectra from FIRS and the SOT/SP. Most pre- 
dominantly noticeable is the effect of seeing on the ob- 
servations. Small scale features with extreme intensity, 
velocity, or polarization signatures are smoothed over in 
the FIRS spectrum, such as the upper line spectrum in 
Figure [2] which is from a position in the light bridge. 
However, the overall behavior of the spectra is similar, 
and non-symmetric profiles seen in the Hinode data are 
faithfully reproduced by FIRS. In the scatter plots the 
overall adherence to the 1:1 line (esp. at low signal levels 
in Stokes I) is encouraging. The flattening of Stokes Q, 
U, and V in Figure [3] indicates that the highest polariza- 
tion signals seen by Hinode are not detected by FIRS, 
these signals are spatially discrete and become smoothed 
over by seeing. The high Q and U values in the Hin- 
ode data where corresponding FIRS values are zero, are 
the result of a bad scan-step on the penumbral/quiet sun 
boundary in the SOT/SP scan. 

3.3. Inversion of Sunspot Magnetic Fields 

We have developed a Milne-Eddington inversion code, 
called the Two- Component Magneto- O ptical (2CM0) 
Inversion Code, based on the formalism of lJefferies et al.l 
|l989) for the inference of the vector magnetic field con- 
figuration from our data. This code models the observed 
Stokes spectra with a magnetic and a non-magnetic 
atmospheric component, fully accounting for magneto- 
optical effects. Voigt and Faraday-Voigt functions are 
used for the synthesis of the Fe I line profiles, while Gaus- 
sian profiles are used for fitting the OH lines. The 2CM0 
inversion code is capable of simultaneously fitting of mul- 



tiple spectral lines with different Lande geff factors, such 
as the line pairs of the Fe I 630 nm and 1565 nm lines, to 
achieve higher ma gnetic field str ength sensitivity in the 
weak-field regime (jSolanki et al.lll992h . 

The magnetic component in 2CM0 is described by ten 
parameters: magnetic filling fraction (/), magnetic field 
strength (5), inclination (7), azimuth (x), source func- 
tion (^o), source function gradient (^i), line center (A), 
doppler width (AA), damping parameter (a), and the ra- 
tio of line to continuum absorption (7^0 )• By assuming 7^0 
has the same behavior between the two lines the num- 
ber of fitted parameters can be decreased. The value 
of 770 is scaled for the atomic lines by a constant factor 
of gife~^^^^^^ where gi is the degeneracy of the lower 
level, / is the oscillator strength, Ei is the energy of the 
lower level for the transition, and T is the temperature 
(we assume an average temperature at the height of line 
formation of 5000 K). The line values log(gf) and Qeff 
for the Fe I lines were taken from Borrero et al. (2003). 

The OH lines are very weak in regions warmer than the 
sunspot umbra, and are therefore susceptible to stray 
light in much the same way as the atomic lines. In 
addition, t he OH lines experience the molecular Zee- 
man effect (Berdvugina fc Solankill2002l ) which produces 
non-negligible Stokes V signal at umbral magnetic field 
strengths. In order to accurately fit the Stokes V pro- 
files during the inversion, and properly account for the 
stray light contamination for the OH lines, they must be 
included in the magnetic component. The OH lines are 
assumed to be infiuenced by the same magnetic field as 
the atomic lines, however the gaussian width and depth 
of the OH lines are left as independent parameters. 

A non-magnetic component is included to account for 
stray light from the instrumental system and unresolved 
mixing of adjacent magnetic and non-magnetic regions. 
It is assumed the non- magnetic fill fraction is 1 — /. A 
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Fig. 4. — Fitted line profiles from the MERLIN (blue) and 2CMO (red) inversions with the original Stokes spectrum from Hinode (black 
points), from the umbra of NOAA 11024. 



stray light profile is generated from averaged quiet sun 
intensity profiles and is allowed to shift in wavelength 
(As), adding one additional free parameter. 

To confirm the validity of the solutions from our inver- 
sion code, we have conducted a comparison of the results 
from 2CM0 and those from the Milne-Eddington gRid 
Linear Inversion Network (MERLIN) from the Commu- 
nity Spectropolarimetric Analysis Center (CSAC) at the 
National Center for Atmospheric Research. MERLIN is 
a new implementation of the ori^nal ME code written 
bv lSkumanich fc LitesI (|1987[ ) (|Lites et al.ll2Q07[ ). MER- 
LIN fits the same set of parameters as our code, and 
although the subtleties of the implementation, initializa- 
tion, and calculation are different, the two codes describe 
the same physics and should produce similar results when 
performed on identical data. 

As an example we have obtained MERLIN results from 
the CSAC inversion client for a 6302 A observation of the 
sunspot in NOAA 11024 taken on 2009 July 7 by the So- 
lar Optical Telescope Spectro-Polarimeter onboard the 
Hinode spacecraft, and preformed a 2CM0 inversion on 
the same dataset. In Figure 2] we show the observed 
Stokes profiles from a single position in the sunspot um- 
bra (black), the fitted profiles from 2CM0 (red) and 
MERLIN (blue), and the data- fit residuals. The syn- 
thetic profiles show slight differences due to the choice 
of line parameters, however the quality of fitting shown 
by the data- fit residuals is comparable. In Figure [5] we 
show the scatter-plot comparison of the magnetic field 
parameters from MERLIN and 2CM0 for all inverted 
positions. The two codes produce largely the same mag- 
netic field strength in the umbra to 70 G (rms). The 
increased scatter in magnetic field strength at low values 
demonstrates the degeneracy with the magnetic filling 
factor in regimes where the line is not completely split. 
However, the scatter in the filling factor is larger than 
we might expect based on this degeneracy alone, and in 



general larger values of the filling factor are produced 
by MERLIN. While the parameter sets are identical be- 
tween the two codes, the implementation is different, in 
particular the treatment of scattered light and the tol- 
erances placed on the Stokes I fit can significantly alter 
the resulting fill factor. The 2CM0 inversion code and 
the results fro m the comparison are described further in 

Q^ldmiD. 

3.4. Treatment of Stray Light 

Stray light is one of the largest sources of er- 
ror in temperature and magnetic field measurements 
in sunspots. Detailed discussions of the effects of 
stray ligh t and its correctio n have bee n carried out 
in .Martinez Fillet fc VazquezI ()1993f ) and iSolanki et al.l 
(|1993f ). In general strav light results in reduced contrast 
of the continuum with respect to the quiet sun, the con- 
tamination of sunspot intensity profiles with a stray quiet 
sun component, and reduced amplitudes of the polarized 
components. The correction of stray light is particularly 
important for the accurate determination of quantities 
which depend on intensity, such as temperature and the 
equivalent width of spectral lines. 

Spatial stray light may arise from both atmospheric 
and instrumental scattered light. Without proper char- 
acterization of the instrumental and atmospheric con- 
tribution to spatial stray light, the fraction due to the 
instrument and telescope cannot be distinguished from 
the unresolved mixing of adjacent magnetic and non- 
magnetic regions. The degeneracy of magnetic field 
strength and its filling fraction in weak field regimes is an 
additional source of ambiguity for the estimation of the 
stray light somponent. Therefore the filling fraction de- 
termined by the inversion can only place an upper limit 
on the spatial stray light. Detailed measurements to de- 
termine the stray light for FIRS have not been carried 
out. However, in observations at the solar limb (shown 
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Fig. 5. — Scatter plot of the magnetic field parameters fitted by MERLIN compared with those from 2CMO from inversion of Hinode 
SOT/SP observations of NOAA 11024 taken on 2009-07-07. Contours are drawn at 52, 84, and 97% levels. 1:1 correspondence is indicated 
by the diagonal black line in each plot. 



1.0 




Fig. 6. Profiles of the 6302 A (top) and 15650 A (bottom) 
continuum intensity at the solar limb as a function of distance 
from solar center. The colored line in each case is the theoretical 
limb darkening function for 0.6 and 1.6 /xm respectively. 

in Figure [6|) the stray fraction at 0.02 R© or about 20" 
above the hmb is approximately 3% and 5% of the disk 
center brightness in the visible and infrared (or about 8% 
of the intensity at the limb). This measurement agrees 
with the filling fraction determined in the umbra of very 
large sunspots, where the photosphere should be almost 
entirely magnetic. Therefore we apply the inverted mag- 
netic filling fraction to correct the spatial stray light in 
the intensity spectrum. 

There is no line- free source of spatial stray light, there- 
fore any spectral stray light must arise only following 
the slit of the instrument. It is possible to estimate the 
spectral stray light by measuring the equivalent width 
of spectral lines and comparing them to a pristine mea- 
surement free of spectral stray light. The presence of a 



possible large spectral stray fraction in FIRS and its ef- 
fect on the measurement of equivalent widths is discussed 
further in M.2\ 

4. INFERENCE OF THE H2 ABUNDANCE 

The direct measurement of H2 at photospheric heights 
is infeasible for the reasons described in ^1.41 Therefore, 
it is necessary to establish the H2 abundance through 
indirect methods. The strength of the OH lines near 
15650 A is strongly related to the temperature at the 
height of continuum formation, although the lines are 
formed higher in the atmosphere by about 100 km. At 
high temperatures the behavior of the lines is adequately 
described by a simple collisional and radiative model 
(iPenn et al.il200 3b). Therefore, they provide us with a 
way to characterize the state of the atmospheric gas in 
addition to the measurement of continuum intensity. We 
establish the equivalent width of OH as a diagnostic for 
the H2 fraction using spectra synthesized from model at- 
mospheres. 

4.1. Model Atmospheres 

In order to characterize the average properties of 
sunspot atmospheres we have examined 1-D atmospheric 
models with solar gravity and abundances, but with a 
range of effective temperatures. We use the Kurucz 
model atmosphere grids (Te//=4000 to 7000 K in 250 
K steps) (http:/ /kurucz. harvard.edu/grids. html), and to 
fully explore the full range of sunspot temperatures, 
cooler models from the updated Phoenix model atmo- 
sphere grids (Te//=2600 to 3900 K in 100 K steps) were 
kindl y provided by Peter Hauschildt (Hauschildt et al.l 
[1993). 

To generate the necessary spectral diagnostics, these 
model atmospheres were supplied to the Rybicky- 
Hummer (RH) code developed by Han Uitenbroek, which 
has been demonstrated in .Uitenbroek (2000a. b). The 
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Fig. 7. Synthetic spectra from the RH code for the FIRS 15650 
A bandpass for each model atmosphere from an effective tempera- 
ture of 2600 K (red) to 7000 K (violet). 
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Fig. 8. — Semi-log plot of fraction of the total particle density 
made up by the H2, CO, OH, and ON molecules at r = 1 in the 
Fe I 15648.5 A line core (black) and continuum (red) as a function 
of the 15650 A continuum temperature and intensity. 

code is able to solve the equations of radiative transfer 
and statistical and chemical equilibrium for atomic and 
molecular species in LTE and non-LTE situations. We 
have used the LTE solutions from the code to generate 
synthetic spectra and obtain the necessary diagnostics 
for our 1-D plane parallel model geometry. In addition 
to calculating the molecular populations at every height 
in the model, we used the RH code to produce a synthetic 
spectrum of the 15650 A range including the Fe I and OH 
lines for each atmospheric model viewed at ten heliocen- 
tric angles (0.0°, 10.2°, 23.4°, 36.2°, 48.5°, 60.0°, 70.3°, 
78.9°, 85.3°, 89.1°). In Figure [3we show the continuum 
normalized spectra from disk center for each atmospheric 
model from an effective temperature of 2,600 K (red) to 
7,000 K (violet). 

The equivalent width of OH and the continuum inten- 
sity (normalized to the quiet sun model Te//=5,750 K) 
were measured directly from the synthetic spectrum for 
each model and heliocentric angle. The corresponding 
molecular abundances were determined at two heights 
in the model atmosphere, r = 1 in the continuum 




4000 5000 6000 7000 

T (T,e=l) [K] 

Fig. 9. — The equivalent width of the synthesized OH 15651.9 
A line plotted against the 15650 A continuum temperature and 
intensity. The H2 density fraction from Figure [8] on a linear scale 
is included for reference and should be read from the right hand 
scale. 
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Fig. 10. — The H2 fraction at the r = 1 height for the the Fe I 
15650 line core (black) and continuum (red) as a function of OH 
equivalent width. 

where the intensity /temperature is measured, and at 
T = 1 in the Fe I line core where the magnetic field 
strength measurement can b e considered to originate 
([Sanchez Almeida et al.l Il996f ). A linear interpolation 
was used to retrieve the parameters for arbitrary helio- 
centric angles and temperatures. 

Figure [8] shows the logarithmic fraction of the total 
particle number of the molecules H2, OH, CO, and CN 
at r = 1 in the 15648.5 A Fe I line core (black lines) and 
continuum (red lines) as a function of the continuum tem- 
perature and intensity at disk center. The OH fraction 
closely tracks the increase of the H2 fraction from about 
6,500 to 3,500 K. Below 3,500 K the OH abundance de- 
creases due to competition with CO for atomic oxygen. 
The other molecular species do not follow the behavior 
of H2 particularly well. This demonstrates the validity of 
OH as a proxy for H2 for atmospheres with temperatures 
above 3,500 K. Figure [9] shows the predicted equivalent 
widths of the OH 15651.9 A line over the range of he- 
liocentric angles spanned by the sample as a function 
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of the 15650 A continuum temperature and intensity at 
disk center (black shaded region). The fractional num- 
ber density of H2 for the Hue core (red dashed line) and 
continuum (black dashed line) are included in linear scale 
to show the good proportionality between OH equivalent 
width and H2 density. Figure [10] shows the direct rela- 
tion between OH equivalent width and the H2 fraction 
at disk center. 

It is necessary to use multiple heliocentric angles to de- 
termine the correct relation between the equivalent width 
of the 15651.9 A OH line and H2. Because the height of 
r = 1 changes as a function of heliocentric angle, the OH 
line strength, intensity ratios, and atmospheric param- 
eters for sunspots at significantly different positions on 
the solar disk cannot be compared directly. To make use 
of observations at obtained over a large range of helio- 
centric angles we must apply a correction based on the 
atmospheric models. For the heliocentric angle of each 
sunspot observation we determine the relation of the in- 
tensity and OH equivalent width at the heliocentric angle 
to the equivalent disk center continuum temperature and 
H2 fraction respectively. In this way data from a wide 
range of angles can be compared directly. 

4.2. Measurement of Eq Won 

We have chosen to use the 15651.9 A OH line which is 
stronger than the blend at 15651.0 A, and less affected 
by the Zeeman split Fe I g=1.7 line than the 15653.4 A 
line. However, the presence of the g=1.7 line still signif- 
icantly affects the inferred equivalent width of the OH 
line at even modest magnetic field strengths. Observa- 
tional systematics in FIRS such as fringing and improper 
correction of the filter profile increase the uncertainty in 
the direct measurement of equivalent width from the in- 
tegrated spectrum. Therefore, in the umbra we calculate 
the equivalent width from the inverted parameters for 
the OH line, and elsewhere in the sunspot the equivalent 
width is calculated directly from the Stokes I spectrum 
because it does not produce poor results when the lines 
are weak. 

Figure [TT] shows the compiled EqWoH vs. continuum 
temperature for the sample. The data have been ar- 
ranged into 50 K bins and the mean and standard devi- 
ation of each bin was calculated. The data and model 
describe qualitatively the same curve which appears to 
be different only by a factor of about 1.55. 

This discrepancy between the empirical and theoretical 
relation of OH equivalent width with temperature must 
come from systematics in the observations or a lack of 
realism in the model atmospheres and the RH code. A 
large amount of spectral stray light could account for the 
consistent underestimation of the OH equivalent width. 
If we assume a model where the stray light level at any 
given position is constant fraction of the quiet sun contin- 
uum {r = I -\- f ^ Iq, e.g. from uniform scattered or par- 
asitic light), then the corrected equivalent width would 
be increased, but the continuum contrast would also be 
increased, and the resulting relation would move up and 
to the left in Figure [Til However, if we assume that the 
stray fraction at any given position is proportional to the 
continuum intensity at that position {I' = / * (1 + /), 
e.g. from overlapping orders), then the equivalent width 
increases by a multiple of the stray fraction, but the con- 
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Fig. 11. Equivalent width of the OH 15651.9 A (red) and Fe I 
15654.8 g=3 (blue) lines compiled from the sample vs. continuum 
temperature at 15650 A are indicated by the points with error bars. 
The theoretical values of the OH and Fe I line equivalent widths 
for heliocentric angles spanning the sample are indicated by the 
solid (disk center) and dashed lines (~ 50°). 

tinuum contrast is preserved. 

If there were a large fraction of spectral stray light it 
should affect the iron lines in the same way. The equiva- 
lent width of the whole Zeeman triplet was calculated for 
the sunspots and the surrounding regions of quiet sun. 
The data was binned and the mean and standard de- 
viation were determined in the same way as for the OH 
line, and is indicated by the blue points in Figure[TTl The 
large scatter in the Fe I measurement makes the possible 
stray light contamination difficult to determine for these 
lines, but it appears that they also underestimate the 
Fe I equivalent width relative to the model. 

We are reasonably confident with the calculation of 
the OH populations and the lines in the synthetic spec- 
tra, however the same cannot be said about the Fe I 
lines. The Fe I lines have non-LTE ionization and we 
have calculated them assuming LTE. The infrared NSO 
sunspot atlas (atlas 5) obtained with the Fourier Trans- 
form Spectrometer (FTS) on the McMath Solar Tele- 
scope at Kitt Peak provides a high-resolution, pristine 
solar spectrum with e xceptionally low spectral stray light 
(| Wallace et al.l [20QTh . From the profiles (available at 
ftp://nsokp.nso.edu/pub/atlas/) we have measured an 
equivalent width of 0.0968 A for the Fe I 15648.5 A g=3 
line which is in very good agreement with the model at 
the quiet sun continuum temperature of 6850 K. In spite 
of this agreement, we can trust neither the measurement 
of Fe I equivalent width, nor the RH model for the same 
quantity beyond the quiet sun, and we must rely on OH. 

The OH and Fe I equivalent widths, adjusted by a fac- 
tor of 1.55 in comparison to the models, are shown in 
Figure [121 If we are to believe this factor, it implies 
a spectral stray light fraction of 35%. This model of 
stray light is not particularly realistic for FIRS, where we 
might expect overlapping orders to fall onto neighboring 
slits. From the favorable comparison of FIRS 6302 A 
with Hinode SOT/SP data, which shows very compara- 
ble contrast in both the continuum and in the lines, we 
can conclude that the possible stray fraction is at least 
limited to infrared. In the future, further verification of 
the RH synthesis of OH is possible using OH lines at visi- 
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Fig. 12. — The same as Figure [TT] but the equivalent widths 
measured from the sample have been scaled by a factor of 1.55. 

ble wavelengths. Additional investigation is necessary to 
confirm the presence of such a large spectral stray light 
fraction in FIRS, but for the moment we apply a uni- 
versal factor of 1.55 to the OH equivalent widths and 
proceed with the H2 inference under this caveat. 

5. SELECTED SURVEY RESULTS 

In this section we provide context for the observa- 
tions using sunspot histories that have been determined 
from SOHO/MDI and SDO/HMI daily intensitygrams, 
and discuss results for each sunspot individually. The 
sunspot observations have been inverted to retrieve the 
magnetic field, the OH equivalent has been measured, 
and model atmospheres have been applied to infer the 
continuum temperature and H2 fraction. 

The results from this analysis for each of our exam- 
ple sunspots are shown in Figures [13] to [19] sorted from 
strongest to weakest by umbral magnetic field strength. 
In panel (a) for each spot observation we show the 
thermal-magnetic relation (square of the magnetic field 
plotted against the continuum brightness temperature), 
umbral (black points) and penumbral (gold points) re- 
gions have been selected based on intensity. In the same 
panel we show the measured OH equivalent width scaled 
by a factor of 1.55 and the OH equivalent width for 
the sunspot heliocentric angle produced by the RH code 
(black line). In panel (b) of the figure we show the in- 
ferred H2 fraction at the height of r = 1 in the Fe I 
15648.5 A line core for a disk center model atmosphere 
plotted against the continuum temperature in blue, the 
theoretical relation between temperature the H2 fraction 
is included in black. Maps of the magnetic field strength 
and continuum intensity are shown in panels (c) and (d) 
respectively. The solid gold line and dotted black line 
in each of the maps show the selected penumbral and 
umbral boundaries respectively. 

NO A A 11131: was a large a-spot that appeared fully 
formed from around the east limb on 1 Dec 2010, 
remained roughly the same size and shape until it 
disappeared around the west limb on 13 Dec. The 
region developed a few transient pores during its 
passage of the disk. The observation in Figure [T3] 
was taken two days before it crossed the central 
meridian. 

NOAA 11131 is the largest spot in the sample with 



TABLE 2 
Summary of Selected Results 



Target Name Max. B [G] Min. T [K] Max. 1^H2/^ Total [%] 



NOAA 11131 


3100 


5100 


2.3 


NOAA 11035 


3000 


5450 


1.5 


NOAA 9429 


2900 


5700 


1.0 


NOAA 11130 


2700 


5650 


0.8 


NOAA 11101 


2600 


5750 


0.9 


NOAA 11046 


2500 


5850 


0.7 


NOAA 11049 


2200 


5950 


0.3 



a maximum umbral field of 3,100 G. The thermal- 
magnetic relation follows a significantly different 
path from the other spots. Large magnetic field 
strengths are reached early in the penumbra fol- 
lowed by a shallow slope through a majority of the 
umbra. An abrupt increase in slope occurs at a 
temperature of 5,300 K. The OH lines are quite 
strong, reaching a value corresponding to 2.3% of 
hydrogen bound into H2. 

NOAA 11035: was a complex /37-region that emerged 
rapidly on 14 Dec 2009, 10° from the central merid- 
ian. Disorganized pore groups coalesced into large 
leading and following polarity sunspots. The obser- 
vation of the leading spot in Figure [M] was taken 
on 17 Dec when the spots were largest. After pro- 
ducing several C-class flares on this day the region 
began to decay, but the spots remained large as 
they traveled around the west limb on 22 Dec. 

The NOAA 11035 observation shows a maximum 
umbral magnetic field strength of 3,000 G in the 
main umbra. The spot shows exceptionally high 
magnetic field strengths in excess of 3,500 G in 
small features located to the upper left and upper 
right of the main umbra (can be clearly identified 
in panel (c)). Interestingly these strong magnetic 
fields appear at the umbra/penumbra boundary of 
the secondary umbrae and do not correspond with 
the coolest temperatures, or strongest OH lines, 
and are probably a manifestation of flare activity. 
Disregarding the scatter in magnetic field caused by 
these discrete magnetic features, the main thermal- 
magnetic relation for this spot does not display a 
sharp upturn, the relation gradually increases in 
slope in the main umbra starting at 5,900 K. The 
H2 fraction inferred from the OH equivalent width 
is in good agreement with the model and reaches a 
maximum value corresponding to a H2 fraction of 
about 1.5%. 

NOAA 9429: was an alpha spot that appeared fully 
formed around the east limb on 13 Apr 2001 and 
disappeared around the west limb on 25 Apr. The 
observation in Figure [15] was taken on 18 Apr just 
before it passed disk center, following this it ap- 
peared to slowly decay. 

This is an early observation taken with the HSG 
prior to the commissioning of FIRS. The NICMOS 
3 detector used for this observation appears to have 
a significantly non-linear response which was not 
measured at the time of observation. However, the 
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Fig. 14. — For the leading polarity spot in NOAA 11035 taken on 2009-12-17 (a) thermal-magnetic relation and OH equivalent width, 
(b) H2 fraction at the height of the Fe I line core inferred from the OH equivalent width vs. the continuum temperature, (c) magnetic field 
map, (d) map of continuum intensity. 



umbral magnetic field strength derived from Stokes 
profile inversion is independent from intensity and 
still a reliable quantity. We have adjusted the in- 
tensity linearly to match the FIRS observations, 
and the intensity-dependent quantities, (temper- 
ature, OH equivalent with, and N(H2)/N(Total)) 
should only be considered in a qualitative way. De- 
spite the uncertainty in the temperature scaling. 



this sunspot shows a striking isothermal increase 
in magnetic field up to strengths of 2,900 G at the 
lowest temperatures. This upturn can be ascribed 
to three strong magnetic cores visible in the center 
of the sunspot umbra seen in panel (c) of Figure 

M 

NOAA 11130: was a rapidly evolving, medium sized 
/3-region with a single lead spot and weaker fol- 
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Fig. 15. — For NOAA 9429 taken on 2001-04-18 (a) thermal-magnetic relation and OH equivalent width, (b) H2 fraction at the height 
of the Fe I line core inferred from the OH equivalent width vs. the continuum temperature, (c) magnetic field map, (d) map of continuum 
intensity. 




Fig. 16. — For NOAA 11130 taken on 2010-12-02 (a) thermal-magnetic relation and OH equivalent width, (b) H2 fraction at the height 
of the Fe I line core inferred from the OH equivalent width vs. the continuum temperature, (c) magnetic field map, (d) map of continuum 
intensity. 



lowing group which emerged near disk center on 
27 Nov 2010. The region reached a maximum size 
on the 30th, the lead spot having developed a full 
penumbra, and following spots with partial penum- 
brae. It slowly decayed over the next 4 days as it 
reached the limb. The observation shown in Figure 
[T6l was taken of the leading polarity sunspot on 2 
Dec during this decay phase when the spot group 



with following polarity had diminished to pores. 

While NOAA 11130 does not have a particularly 
strong magnetic field (2,600 G) or a large H2 frac- 
tion, this sunspot also shows a clear case of isother- 
mal intensification of the umbral magnetic field. 
The magnetic field turns sharply upward, increas- 
ing by approximately 500 G starting at 5,800 K. 
The inferred H2 fraction for this sunspot in particu- 
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Fig. 17. — For NOAA 11101 taken on 2010-09-02 (a) thermal-magnetic relation and OH equivalent width, (b) H2 fraction at the height 
of the Fe I line core inferred from the OH equivalent width vs. the continuum temperature, (c) magnetic field map, (d) map of continuum 
intensity. 
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Fig. 18. — For NOAA 11046 taken on 2010-02-13 (a) thermal-magnetic relation and OH equivalent width, (b) H2 fraction at the height 
of the Fe I line core inferred from the OH equivalent width vs. the continuum temperature, (c) magnetic field map, (d) map of continuum 
intensity. 



lar is in rather poor agreement with the theoretical 
prediction, with a maximum of 0.8% H2. 

NOAA 11101: appeared around the east solar limb as 
a fully formed a-spot with a penumbra on 25 Aug 
2010 and did not experience significant evolution 
as it transited the disk and disappeared around the 
west limb on 5 Sep. The observation shown in Fig- 
urefTTlwas taken just after it had passed disk center 



on 2 Sep. 

This sunspot shows a large scatter in the thermal- 
magnetic relation from an extended magnetic 
structure which spans the umbra-penumbra bound- 
ary on the lower right hand side of the sunspot um- 
bra. However, it shows a higher density of points 
at 5,800 K rising isothermally out of the main 
B^-T track. Agreement with the theoretical H2- 
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Fig. 19. — For NOAA 11049 taken on 2010-02-19 (a) thermal-magnetic relation and OH equivalent width, (b) H2 fraction at the height 
of the Fe I line core inferred from the OH equivalent width vs. the continuum temperature, (c) magnetic field map, (d) map of continuum 
intensity. 



temperature relation is excellent and we can infer 
a maximum H2 fraction of 0.9%. 

NOAA 11046: appeared around the east limb as a 
small /3-region on 7 Feb 2010. It became more de- 
veloped around 11 or 12 Feb with a larger follow- 
ing polarity group and a smaller leading polarity 
sunspot. The observation in Figure [T8l is of the fol- 
lowing spot group on the 13th as the region began 
to decay over the next 5 days. The spots seemed 
to have disappeared almost entirely as it reached 
the limb on 18 Feb. 

Although this sunspot group is composed of three 
umbrae with partial penumbrae, it displays a cohe- 
sive thermal-magnetic relation which shows a clear 
increase in slope in the umbra, at about 6,100 K. 
A maximum umbral field strength of 2,500 G is 
reached in the right-most umbra. From the OH 
measurement we can infer a maximum H2 fraction 
of just under 0.7%. 

NOAA 11049: was a small /3-region that emerged on 
16 Feb 2010 about 15° east of the central meridian. 
Although the region continued to evolve, the lead 
and following pores remained roughly equivalent 
in size until 21 Feb when the lead pore became 
stronger and the following region diminished before 
the region disappeared around the west limb on 24 
Feb. The observation shown in Figure [19] of the 
following polarity pore group was taken on 19 Feb 
during this relatively quiescent evolution stage. 

This pore group shows a very linear thermal- 
magnetic relation and low magnetic field strengths, 
a 2,200 G maximum is reached in the left-most 
pore. The OH lines are small or undetected, corre- 
sponding to 0.3% H2 fraction at most. 



For the sample as a whole, the overall agreement of 
the scaled OH equivalent width with the theoretical OH 
curve, especially in the larger sunspot s, strengthens our 
confidence in the results of the radiative transfer and 
chemical equilibrium code. The maximum umbral field 
strength occurs at the same location as the maximum OH 
equivalent width (except for the fiaring region, NOAA 
11035) as we would expect if OH is strongly linked to 
atmospheric temperature, and temperature is controlled 
directly by magnetic suppression of the transport of con- 
vective energy. 

Among the sample, the medium-sized spots NOAA 
9429 and 11130 show clear isothermal increases in the 
magnetic field strength which occur around T=5,800 K. 
NOAA 11101 shows a larger scatter of the points in the 
vs. T curve, with a higher density track similar to 
those of NOAA 9429 and 11130, it is suggestive of the 
existence of isothermal magnetic field intensification pro- 
cess. For the larger spots, NOAA 11131 and 11035, and 
for the smaller fragmented sunspot group NOAA 11046, 
a clear case for isothermal magnetic field intensification 
cannot be made. However, it should be noted that all 
of the larger sunspots show some kind of increase in 
the slope of the thermal-magnetic relation starting below 
6000 K where the H2 fraction begins to increase rapidly. 
The pores in NOAA 11049 show strictly linear behavior, 
indicating that for the case of small pores, the magnetic 
field is too weak and the temperature too high to have a 
significant H2 fraction, which is confirmed by the unde- 
tected OH lines. 

A very intriguing result of this survey is the vs. T 
curve of the large spot of NOAA 11131 that is distinctly 
different from those of the rest of the sample. It follows 
an elevated track in B^ vs. T without the signature of 
isothermal magnetic field intensification at T = 5,800 K. 
Similar results are found in the other large sunspots in 
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the sample. 

6. CONCLUSIONS AND DISCUSSION 

Our survey and analysis provide observational evi- 
dence that significant H2 molecule formation is present 
in sunspots that are able to maintain maximum fields of 
greater than 2,500 G. Measurements of the OH equiva- 
lent width seen in sunspot umbrae are qualitatively con- 
sistent with the predictions from spectra synthesized by 
radiative transfer models. We infer a molecular gas frac- 
tion of a few percent H2 in the largest sunspots. The 
formation of this small faction appears to alter the equi- 
librium of pressures in the sunspot isothermally, resulting 
in an increase of the slope of the thermal-magnetic rela- 
tion at temperatures lower than a 15650 A continuum 
brightness temperature of 6,000 K where the H2 fraction 
begins to rapidly increase with temperature. We suggest 
that the formation of H2 molecules in the sunspot umbra 
causes a rapid intensification of the magnetic field with- 
out a significant decrease in temperature which would 
explain the increase in slope of the thermal-magnetic re- 
lation. 

We hypothesize that H2 plays an important role in the 
formation and evolution of sunspots. During the initial 
stage of sunspot emergence and cooling, the formation 
of H2 may trigger a temporary "runaway" magnetic field 
intensification process. As magnetic fiux emerges and 
strengthens, the sunspot atmosphere cools due to the 
suppression of convective heating by the magnetic field. 
When sufficiently low temperatures are reached H2 be- 
gins to form in substantial numbers in the coolest parts of 
the umbra. As free hydrogen atoms combine to form H2, 
the total particle number density is reduced. The dissoci- 
ation energy released into the atmosphere is rapidly dis- 
sipated by radiative cooling due to the low opacity of the 
photosphere, reducing the total kinetic pressure without 
a corresponding reduction in temperature. The decrease 
in gas pressure causes this region to shrink in size, and 
due to the high electrical conductivity of the atmosphere 
the magnetic fields are compressed with the plasma (the 
"frozen-in field" effect). The resulting higher magnetic 
field density further inhibits the convective heating of 
the sunspot atmosphere, which leads to further cooling. 
This "runaway" magnetic field intensification process is 
most likely a temporary phenomenon which is arrested 
long before all the hydrogen atoms have condensed into 
molecular form. At some point the transport of energy 
by convection will be effectively quenched and increases 
in magnetic field will no longer result in decreases in 
temperature. The transfer of radiative energy from sur- 
rounding hotter regions would also keep the umbra from 
becoming excessively cool. Therefore further H2 forma- 
tion would be halted. 

While the formation of H2 may initiate a more rapid in- 
tensification of the sunspot magnetic field during sunspot 
emergence, we speculate that during the decay phase in 
a sunspot which has already formed a substantial H2 
population, the highest concentrations of molecular gas 
would tend to maintain the magnetic field against decay 
and extend the lifetime of the sunspot. As the magnetic 
field in a sunspot weakens, regions of the umbra once 
cool become warm and H2 dissociates back into atomic 
hydrogen. The more rapid increase in pressure in warmer 
regions of the umbra would compress the remaining cool 



regions, concentrating the magnetic field and maintain- 
ing the cool interior against convective heating. 

The formation of H2 would speed up the process of 
sunspot emergence, and the dissociation of H2 would slow 
down their disappearance. While the effects of the forma- 
tion and destruction of molecules would produce similar 
signatures in the thermal-magnetic relation, we would 
expect to see more cases of sunspots in the decay phase 
due to observational bias. There is evidence that the in- 
tensification of the magnetic field occurs in discrete cores, 
such as can been seen in NO A A 9429, which is consis- 
tent with this speculative picture of molecule formation 
during the growth and decay of sunspots. If this is the 
case the molecular fraction would be significantly under- 
estimated in sunspots due to the effect of filling factor, 
and may exist in quantities of 5% in unresolved features, 
consistent with the coolest mode ls of the umbra such as 
those presented in Malt by et al.l (0^986). 

It is possible that two other effects contribute 
to the non-linearity of t he ther mal-ma gnetic rela- 
tion. In previous studies (iMartin ez Pillet fc VazquezI 
[T993I : ISolanki et all [1991 iMathew et al. 200^ the non- 
linearity of the thermal-magnetic relation was inter- 
preted as a radiative transfer effect, i.e. the Wilson 
Depression effect in sunspots. Cooler atmospheres in 
a sunspot are more optically thin, therefore the mag- 
netic field and continuum measurements originate from 
a greater geometrical depth. In seeing deeper into the 
atmosphere we are able to see relatively hotter regions, 
therefore the observed temperature would seem to de- 
crease less rapidly (relative to radius or B) than in a sin- 
gle geometrical layer. This effect would therefore tend to 
increase the slope of the thermal-magnetic relation, how- 
ever the temperature should still decrease as the mag- 
netic field strength increases. 

Through all of this work we have also neglected the 
curvature force. Increased contributions to the horizontal 
support from the curvature force in outlying regions may 
cause the magnetic pressure in the sunspot core to seem 
boosted, contributing to the non-linearity of the thermal- 
magnetic relation. 

Detailed modeling efforts are necessary to determine 
the vahdity of the proposed scenarios for II2 forma- 
tion and destruction during the emergence and decay of 
sunspots, and contribution of the radiative transfer ef- 
fect and the curvature force to the non-linearity of the 
thermal-magnetic relation. We intend to more carefully 
consider the effects of the curvature force and the Wilson 
Depression and further investigate the problem of H2 for- 
mation using the simultaneous dual-height observations 
obtained with the 6302 and 15650 A channels of FIRS to 
perform a detai led co mparison with recent MHD sunspot 
models from iRempell (j2010). 

Particularly for the cases of isothermal intensification 
of the magnetic field in NOAA 9429 and 11130, it is 
unlikely that such a sharp upturn in the slope can be 
explained simply through radiative phenomenon or the 
curvature force. The formation of H2 appears to be the 
most likely cause of the sharp increase in slope of the 
thermal-magnetic relation in sunspots. In addition to 
this magnetic intensification process, the formation of 
molecules increases the heat capacity of the sunspot at- 
mosphere. Due to the additional degrees of freedom of 
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the H2 molecule, the formation of a H2 fraction of 10% 
would ideally raise the heat capacity of the gas by 13% 
over an equivalent number density of atomic gas. This 
non-thermal reservoir for energy may have an important 
effect on the local radiative output of the Sun. Conse- 
quently, we suggest that modeling of the MHS equilib- 
rium condition of sunspots in the form of Equation [1] 
must include a multiple-component atmospheric model 
with the proper equation of state to account for the al- 
tered thermodynamics of the sunspot atmosphere due to 
the formation of H2 . 

Our sample was obtained sporadically through the 
minimum phase of solar cycle 23 and does not contain 
very large sunspots with high magnetic field strengths. 
An intriguing and unexpected finding is the distinctly 
different behavior of the vs. T curve of the largest 
sunspots in the survey (e.g., NOAA 11131). While we 
cannot provide a definitive explanation at this point, we 
point out that Equation 1 describes only the MHS equi- 
librium condition for sunspots, and sunspots cannot be 
in perfect MHS equilibrium all the time (in such case 
sunspots would be static without possibility of evolu- 
tion). Therefore we should not expect that the vs. 
T curves to represent the equilibrium state, or that they 
follow the same track. We suspect that differences in the 
observed behavior of the B^ vs. T curves are the manifes- 
tation of changing magnetic and thermal environment of 
sunspots at different stages of their evolution. A contin- 
ued observational effort following sunspots through their 
Hfecycle should provide the necessary data to address 
this issue as solar cycle 24 enters its maximum phase 
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